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SUMMARY 


This note presents a study of the current methods of computing the 
Doppler observable and, specifically, a range-difference technique for 
computing integrated Doppler count. The integrated Doppler count can be 
directly computed from the difference in lengths of the paths traversed 
by the transmitted signals at the beginning and end of the counting 
interval. When determining rajectories in the vicinity of the earth- 
moon system, this range-difference technique has certain advantages over 
an existing technique called the range-rate expansion method. There is 
no truncation error for the counting intervals presently planned for Apollo 
missions, and the formulation and resulting implementation are much shorter 
and less complex. Thus, this range-difference technique for computing the 
Doppler observable waa chosen for use in the Apollo orbit determination 
program. 

The derivation of the Doppler equation presented in this report employs 
Newtonian methods. In the computation of the position vectors of the 
vehicle and stations, the finite value of the velocity of light in the 
medium is taken into account. This approach is valid since, for two-way 
Doppler, the effects of general and special relativity on the observable 
in the earth-moon region are less than the known noise on the measurement; 
for three-way Doppler, these effects are significant only for widely sepa- 
rated stations and produce a small bias which may be solved for. 

A numerical study was performed in which the total number of Doppler 
counts over the counting interval t was computed using the range-rate 
expansion method. The truncation error was reduced by dividing the inter- 
val into smaller subintervals , computing the number of Doppler counts over 
each subinterval, and summing the results. This was taken as the standard 
value. From the same trajectory, at the same times, the Doppler count was 
computed using the range-difference method. The difference in the two 
values of Doppler count was then recorded. 
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Three cases were studied: earth orbit with a maximum elevation of 

80° , eariy tranalunar, and lunar orbit. In the earth orbit case the maximum 
difference between the two methods was 6.0 * 10 5 cps or 1.1* * 10 5 fps_in 
range-rate. _ In the translunar case the maximum difference was 2.0 * 10 3 cps 
or 0.5 x 10~ 3 fps. In the lunar orbit case both two-way and_three-way 
Doppler were studied and the maximum difference was 3.0 x 10 5 cps or 
0.7 * 1C” 5 fps. These values are orders of magnitude less than the instru- 
ment noise which for 6-second counting intervals is approximately 0.2 cps 
cr 0.0U fps in range-rate. 

These results show that, when the truncation error of the range-rate 
method has been removed, the two methods agree to well within measurement 
capabilities . 

INTRODUCTION 

Unified S-band Doppler count is the most fundamental and important 
ground-based observation used during Apollo missions for locating the 
spacecraft. The use of Doppler observable can be briefly described as 
follows. The Doppler count is recorded by the tracking stations. The 
real-time orbit determination program uses its current estimate of the 
vehicle’s position and velocity to compute a corresponding value of Doppler 
count. The difference between the observed and computed values of Doppler 
count, called the residual, is compute* and used by the orbit determination 
program to obtain an improved estimate of the vehicle’s position and 
velocity. 

The basic, continuous count Doppler observable is illustrated in 
figure 1. Two- and three-way Doppler geometry is shown. The transmitting 



Figure 1.- Doppler geometry. 
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station transmits the frequency v,p. The vehicle receives the signal, 

multiplies it by a constant, , and retransmits. The receiving station, 

which is the same as the transmitting station in the case of two-way Doppler 
and is a different station in the case of three-way Doppler, records the 
accumulated Doppler count. 

The contents of this counter are recorded at precise intervals of 
time (usually either 6-second or 60-second counting intervals), and are 
transmitted to the real-time orbit determination program. After receiving 
this information, the difference between successive values of count is 
formed, and the result is normalized by dividing by the counting interval. 
This, then, is the expression that must be computed by the orbit deter- 
mination program so that a residual can be formed, and an improved estimate 
of the vehicle trajectory obtained. 


Range-Rate Expansion Method 

To compute the expression, the integrated Doppler count over the 
counting interval, t, must be computed. This can be accomplished by 
expanding the frequency shift in a power series of t, and integrating 
term by term, evaDuating the coefficients in the middle of the counting 
interval (ref. 1, 2, and 3). In the equations in references 1 and 2, only 
the first four terms of the expansion (i.e., the i° and t 2 terms) are 
carried (the t 1 and t 3 terms are zero). 

It is pointed out in reference k that this power series expansion is 
adequate for presently planned x intervals during the late trans lunar and 
lunar phases of Apollo type missions. However, it is also pointed out 
that in certain worst geometrical situations occurring in earth orbit and 
immediately after translunar injection, the errors resulting from the 
truncation in the t series are significant. 

This truncation error is illustrated in figure 2. In this case, an 
earth orbit pass with a maximum elevation of 87° is shown. The truncation 
error, expressed both in cycles per second and feet per second, is plotted 
versus elevation angle. The tracking station instrument errors are approxi- 
mately 0.2 cps for a 6-second counting interval, and it can be seen that 
between the 60° elevation points, the truncation error is larger than the 
instrument error, 
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Figure 2.- Truncation error. 


One method of avoiding the truncation error in these cases is to 
break the x interval into several smaller intervals. For example, a 
6-second t interval could be broken into six, 1-second x intervals. The 
expansion could then be applied to each of the smaller intervals, and the 
results summed to obtain the result for the entire x interval. This, how 
ever, introduces mechanization problems into the orbit determination 
program and adds considerably to the computation time. 


Range-Difference Method 

Another method of computing the Doppler observable is discussed in 
references 5 through 9. This range-difference technique computes the 
integrated Doppler count directly from the difference in lengths of the 
paths traversed by the transmitted signals at the beginning and end of 
the counting interval. 

There are two main advantages of the range-difference method over the 
previously discussed technique. First, since no expansion in powers of 
x is made, no truncation errors proportional to the length of x exist. 
Secondly, the formulation and resulting implementation are much shorter 
and less complex than is the case with the approach which uses an expansion 
in x. 
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The main di sadvantage of this approach is that, in lunar orbit, 
several significant digits are lost when a range-difference computation 
is made. This is due to the large range values, of the order U x 1C 10 cm, 
which exist in lunar orbit. However, when the computations are performed 
in double precision this truncation effect in the range-difference equation 
is not significant. 

A derivation of the equations for computing integrated Doppler count 
using the range-difference technique is discussed in detail in reference 9. 
A more direct approach is presented in this report. In both developments 
of the Doppler equation, Newtonian methods are employed and the finite 
value of the velocity of light in the medium is taken into account when 
computing the position vectors of the vehicle and stations. 

This approach is valid since for two-way Doppler, the effects of 
general and special relativity on the observable in the earth-moon region, 
and for vehicle velocities involved, are less than the known noise on the 
measurement, and for three-way Doppler these effects are significant only 
for widely separated stations and produce a small bias which may be solved 
for. 


In view of the advantages of the range-difference technique over the 
method which uses the expansion in t and since the Apollo orbit determina- 
tion program is in double precision, the range-difference technique for 
computing the Doppler observable was chosen for U3e in the Apollo orbit 
determination programs . 

The main text of this note presents the derivation of the Doppler 
equation. The light-time delay iterative technique, which is used to 
compute the position vectors of the vehicle and the transmitting station, 
is then discussed. Finally, the partial derivatives of the Doppler observ- 
able with respect tc the 3tate and biases are .'resented. 



SYMBOLS 

unknown Doppler bias which may be solved for 
speed of light 

integrated Doppler count observable 

accumulated Doppler count from time t to time t_, 

o R 

accumulated Doppler count from time t to time t -t 

o R 
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s. 



Btate rector of the receiving station in earth-centered 
inertial coordinates 



state • ;tor of the transmitting station in earth-centered 
inertial coordinates 




u) M , U)’ 




vehicle state vector in earth-centered inertial coordinates 
obtained from the vehicle ephemeris 


the last time at which the Doppler counter at the receiving 
station was initialize! 

time at which the signal was received at the receiving 
station 

time at which the signal left the transmitting station 

time at which the nignal arrived at, and was transmitted 
by the vehicle 

propagation time for signal which reaches receiving 
station at time t 

instantaneous phase 

instantaneous received frequency 

instantaneous ti.uismitted frequency 

range value from the transmxtter to the vehicle 

range value from the receiver to the vehicle 

Doppler count sampling interval 

the Doppler bias which will result, in the case of three- 
way Doppler, when the standard frequencies of the two 
stations are set to different values 

a large positive constant frequency, added to the incoming 
frequency at the receiving station, which is designed to 
keep the Doppler count positive over the sampling interval 

a frequency multiplying factor by which the transponder 
in the spacecraft alters the frequency received by the 
spacecraft. 
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DERIVATION OF DOPPLER EQUATION 

A derivation of the Doppler equation is presented in reference 9. 
However, the development presented below is more straight forward and more 
closely describes the functions of the Doppler count measuring equipment 
at the receiving station. 


The analytical expression which describes the accumulated Doppler 
count, D , measured at the receiving station is 


V T " V FEC) dt 


The difference between successive values of count, D. - D , 

l t i-l 

is formed in the orbit determination program. This is then normalized by 
dividing by the counting interval 1 * ^ ^ ^ Thus 


D - D 

t t -T 

R R 


w 3 ♦ b + w' + i f R o u L - VrJrec dt, (2) 

T ' 


where the values v T and are instantaneous values of frequency and 

both are expressed at the same instant of time, measured at the receiving 
station, and the counting interval, x, begins at time t R -x and ends at 

time t R . 

The frequency can be expressed in terms of phase by the expression 


Equations (2) and (3) can be combined to form the following expression 


D = + b + 




d< J> T d4> REC 


dt /REC 
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Equation (U) can be written as 


D ■ ♦ b ♦ 


/ ,t R 

♦ rf *( 

•'t^-T X 


" *rec)rec* 


An exact evaluation of the integral in equation ( 5 ) can be obtained. 
This results in equation (6) which expresses the integrated Doppler count 
in terms of the instantaneous phase difference as measured at the receiving 
station at the beginning and ending times of the counting interval. 

D ■ 103 + b + u' + — j^ T - ♦REC^t R * ( l> T ” ^REC^tp-tJ . (6) 

The instantaneous phase difference must be evaluated in terms of 
tracker, vehicle geometry. 

For a particular receiver time, t, it has taken a finite amount of 
time for the signal to have traveled from the transmitter to the receiver. 
It is assumed that phase propagates, unaltered, through the medium. 

At time, t, the receiver phase will, therefore, differ from the trans- 
mitter phase by the amount that the transmitter phase has changed during 
the propagation time. 

In equation 7 , it can be seen that the amount that the transmitter 
phase has changed is given by the transmitter frequency, v T , multiplied 

by the propagation time At^ , 


♦t “ *REC 


:)t = v 


T At t 


The propagation time is equal to the total distance, p^, t; aveled ^y 

the signal which reaches the receiver at time t, divided by the velocity 
of light 


At. « — 
t c 
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The total distance traveled by the signal is the sum of the uplink 
distance, p^, and the downlink distant, p^, 


t ■ ( p i * 0 2 )t 


Combining these expressions, one obtains the instantaneous phase 
difference in terms of the transmitted frequency and the up and downlink 
distances that were traveled by the signal which arrived at the receiver 
at time t. 


P, ♦ P, 


rT " *REC 


\ . v P ‘'2/t 
:lt t c 


Substituting equation (8) into equation (6), the final result can be 
obtained which expresses the Doppler observable in terms of the paths 
traveled by the signal at the beginning and ending times of the counting 
interval , 




K' 1 * ' 2 )*r ■ ( r * * ° 2 k->] 


. (9) 


The method used to compute the distance traveled by the signal is 
presented in the next section. 


LIGHT-TIME DELAY 


Figure 3 shows a signal which leaves the transmitter at time t^,, 

travels a distance p_ , and reaches the vehicle at time t . The vehicle 

1 v 



Figure 3.- Light time-delay geometry. 
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p 2 * iyv - yy 


* i^ v (^ v ) ~ ^(t^) i 


(lla) 

(lib) 


For the first iteration the position vectors R v , Rp, and R,p are 
computed for time tp. These values are used to obtain and p^, and t^ 

and t^, are then computed from equation 1C. 

These values of t and t_ are used to obtain better estimates of R 

v T v 

and R^, and, hence, better estimates of p^ and p^. Thjs iterative process 

continues until both t v and t T are within specified tolerances. The vector 

R v is obtained from the vehicle ephemeris using Lagrange interpolation. A. 

detailed discussion of this iterative procedure is presented in 
reference 10. 


This method is used to compute the total distance traveled by the 
signals which reach the receiving station at the beginning and ending 
times of the counting interval. These values are then used in equation (9) 
to obtain the integrated Doppler count observable. 
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DOPPLER PARTIALS 

To compute the Doppler partials the assumption is made that 


(»! * - (»l * »*) V ' 




( 12 ) 


Using the light-time delay technique, the vehicle ephemeris , and the 
position vector of the receiving station at t R - ^ as input, the following 
quantities are computed: 

’!• p 2* - r - i) • K^* - r - i) • 

To simplify the terminology, the time subscripts will be omitted. 

Using the approximation made in equation (12), equation (9) becomes 


D = u> 3 ♦ b + 


or 


D * u>^ + b + a)* + 


+ -^r [k ( p i + p s)] 
v» (*> i , dp g\ 

c l dt dt / 


(13) 


Now = |R v - R,p| , and for any vector A = aA, we have the relationship 


. dA da _ 
A * dt " a dt • 


dP 

dt 


1 (R v - R T ) • ( K - V R v * ( R v - V - R T • 


(R v - V 


(lUa) 
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and, similarly. 


do 2 (K v - R r ) • (R v - R r ) R v 

• < R V - 

dt P 2 


From equation (13) we have 


3D ( Vt 

3 / dp l 

3R = c 

3R l dt 

V 

v \ 


(lUb) 


(15) 


From equation (lit), we have 

3 / dp l\ "y- K , - l(R v- V • V * P 1 
aR vV dt /" P 1 *? 3R v ' 


Now 



and, therefore, 


3P X R 



[(% - ty * < R y - R T ) j 

r (R v • R v - 2R v • + "T • V 

- «T 

a 
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Therefore, we have 



and using equation (lU), this becomes 



Similarly, it can be shown that 



1 ^2 
o? dt 




) 


Substituting these equations into equation (15) we have 





1 dp 2 
1 


OX dt 




From equation ( lU ) we see that 




P 


1 


and 


Ih 


Therefore, we have the result that 



<Vt 




It can be shown that these results and the approximations involved 
are the same as those given in reference 3. 

From equation (9), it is seen that the partial of the observation 
with respect to the bias is given by 


3D 

3b 
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NUMERICAL RESULTS 


The procedure used to numerically verify the range-difference method 
is shown below. 




N* 

P 


3 


i — i — i — i — t 
i — 




Range-rate method 


Range-difference method 


The total number of Doppler counts over the counting interval t , 
N*=N« + N* +N» +...+N* was computed using the range-rate 

p p x p 2 p 3 p n * 

expansion method. The truncation error was reduced by dividing the 
interval into smaller subintervals , computing the number of Doppler counts 

over each subinterval (N« , N* N* ) , and summing the results. 

P 1 p 2 p n 
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This was taken as the standard value. From the same trajectory at the 
same times, the Doppler count was computed using the range-difference 

method (N a ). The difference in the two values of Doppler, AN = N* - N 

Ap P Ap 

was then recorded. 

Three trajectories were studied: earth orbit with a maximum elevation 

of 80°, early translunar, and lunar oibit. In the lunar orbit case both 
two-way and three-way Doppler were studied. 

The maximum values of AN for these cases are presented in the follow- 
ing table. 


Case 

Difference, cps (fps) 

mmm 

n 

Eartn orbit 

6.0 X 10" 5 (l.U X 10" 5 ) 


60 

Translunar orbit 

2.0 x 10 -3 (0.5 * 10" 3 ) 


6 

Lunar orbit 
(2-way and 3-way) 

3.0 x io -5 (0.7 x 10 -5 ) 

60.0 

60 


In all cases the values of AN are orders of magnitude less than the 
instrument noise, which for 6-second counting intervals, is approximately 
0.2 cps or .0U fps in range-rate. 

The truncation error in the range-rate Doppler equation was also 
studied. This was done by comparing the Doppler count obtained by summing 
six consecutive 1-second count values, with the count obtained by using 
a counting interval of 6 seconds. In both cases, Doppler was computed 
using the range-rate equation. An earth orbit trajectory, with a maximum 
elevation of 87°, was used. The truncation error, illustrated in figure 2, 
is expressed both in cycles per second and feet per second and is plotted 
versus elevation angle. From figure 2, it can be seen that, the trunca- 
tion error, between the 60° elevation points, is larger than the instrument 
error (0.2 cps ) . 


CONCLUSIONS 

Two conclusions can be drawn from the numerical results. First, a 
significant truncation error exists in the range-rate Doppler equation for 
certain earth orbit geometry situations when presently planned Apollo 
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mission counting intervals are used. Second, when the truncation error 
of the range-rate method has been removed, the range difference Doppler 
agrees with the adjusted range-rate value to well within measurement 
capabilities . 

The advantages of the range-difference method are 

1. No truncation error exists. 

2. Formulation and resulting implementation are less complex than the 
range rate expansion method. 

3. This method gives more physical insight into the Doppler ob- 
servable and the various phenomena which affect it, e.g., atmospheric 
refraction, station location erro" j, and time tagging errors. 

The limitations of the range-difference method are 

1. Computations must be performed in double precision for lunar 
distances . 


2. Double precision would be marginal for deep space missions. 

In view of the advantages of the range-difference method, it was 
chosen for use in the Apollo real-time orbit determination program. 
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